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ABSTRACT 


A viscoelastic  model  is  developed  to  describe  the  behavior  of  a polyurethane  rubber  provided  to 
NIST  by  NSWC,  Carderock.  The  model  is  based  on  concepts  of  time-dependent  rubber  elasticity 
and  single  integral  models  of  viscoelasticity.  The  strain  energy  density  function  is  assumed  to  be 
of  a Valanis-Landel  form,  while  its  time-dependence  is  analyzed  using  the  “perfect  elastic  fluid” 
model  of  Bernstein,  Kearsley  and  Zapas.  In  addition,  time-strain  separability  and  time- 
temperature  superposition  are  assumed  to  apply  to  the  polyurethane  in  the  temperature  range  of 
interest.  Single  step  stress  relaxation  experiments  in  uniaxial  extension  are  used  to  determine  the 
model  parameters.  From  these,  multiple  step  stress  relaxation  and  creep  histories  are  applied  to 
the  material  and  it  is  found  that  the  model  can  describe  the  nonlinear  time-dependent  viscoelastic 
response  of  the  polyurethane  resin  to  these  cyclic  histories  to  within  approximately  ±10  % of  the 
load  or  deformation  value  applied. 
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I.  INTRODUCTION 


The  isothermal,  nonlinear  viscoelastic  response  of  polymeric  materials  above  the  glass 
transition  temperature  has  been  relatively  well  described  using  the  perfect  elastic  fluid  model  of 
Bernstein,  Kearsley  and  Zapas  (BKZ)[1,2].  This  is  particularly  true  when  one  considers  loading 
histories  in  which  the  deformations  or  loads  of  interest  are  monotone  increasing  [3].  In  addition 
to  giving  a reasonable  representation  of  the  nonlinear  viscoelastic  response  of  polymer  melts,  the 
BKZ  model  has  some  advantages  of  simplicity.  The  response  functions  are  determined  from 
single  step  stress  relaxation  experiments  [4]  and  the  formalism  looks  much  like  a “time- 
dependent”  elasticity.  One  drawback  to  the  BKZ  model,  however,  has  been  the  fact  that  it  is 
formulated  in  terms  of  the  invariants  of  the  (relative)  deformation  tensor  with  the  result  that  the 
material  parameters  for  each  geometry  of  deformation  must  be  obtained  experimentally  from 
stress  relaxation  experiments  in  that  geometry.  The  result  is  a need  for  a large  number  of 
experiments  to  make  a full  three-dimensional  model  of  material  behavior. 

Recently,  however,  there  has  been  work  carried  out  in  this  laboratory  [5,6]  in  which  a 
hybrid  Valanis-Landel  (VL)[7]-BKZ  model  was  developed  and  used  to  predict  the  response  of  a 
poly(vinyl  chloride)  melt  to  multiaxial,  constant  rate  of  deformation  experiments  from  single  step 
stress  relaxation  experiments  in  uniaxial  extension.  It  is  the  goal  here  to  further  develop  this 
approach  with  the  aim  of  describing  the  long  term  creep  response  of  polyurethane  materials  in 
use  by  the  U.S.  Navy  from  single  step  stress  relaxation  experiments  in  uniaxial  extension.  The 
work  described  in  this  progress  report  is  that  performed  during  the  period  from  September,  1996 
through  April,  1998  with  the  purpose  of  evaluating  the  utility  of  the  hybrid  VL-BKZ  model  as  a 
tool  for  estimating  the  long  term  creep  performance  of  polyurethane  rubber  samples  provided  by 
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the  NSWC,  Carderock  to  NIST.  Subsequent  work  will  consider  and  model  moisture  effects  on 


the  resin  behavior. 

II.  METHODS  OF  ANALYSIS 

Large  Deformations 

The  states  of  strain  in  large  deformations  are  commonly  described  by  the  principal 
extension  ratios,  AJt  A2,  As  or  by  three  strain  invariants  whose  values  are  independent  of  the 
coordinate  system  [8], 

3 3 


(1) 


The  polyurethane  material  studied  here  is  assumed  to  be  incompressible,  which 
introduces  the  further  relation: 


I,  A.  i 1 

3 1 2 j 


(2) 


The  Valanis-Landel  Strain  Energy  Function  for  an  Elastic  Material 


For  an  incompressible  material,  Valanis  and  Landel  [7]  postulated  the  existence  of  a 
strain  energy  function  W(I,,  I2, 1J  = W(A},  A2,  AJ  that  is  separable  in  the  principal  stretches  A.,,  k2. 


W(XVX2,X3)  = w(A,j)  + w(X2)  + w(Z3)  + a\rt{XxX2X^) 


(3) 


where  w(A)  will  be  referred  to  as  the  VL  function  and  a is  an  arbitrary  constant  that  can  be  set  to 
0 for  an  incompressible  material. 

The  principal  stress  differences  can  be  written  in  terms  of  the  derivative  w ’(A)  of  the  VL 
function: 

dX  (4) 

G1  - oy  = X.w'iX)  - Xjw'(Xp 

For  a uniaxial  deformation,  the  above  equation  can  be  reduced  to: 

Xl=X2 =—  ( incompressibility ) 

"A 

(5) 

with  A a = an  - : A a(A)  = Xw'CX)  - — w/(— ) 

v/i  yi 

w ’(A)  can  be  obtained  from  a combination  of  data  from  tension  and  compression 
measurements  using  the  approach  of  Kearsley  and  Zapas  [9]  using  a recursive  set  of  equations. 
By  forming  the  sequence  of  equations  of  the  form  of  equation  5 , replacing  X by  Aj. , and  by 
summing  them,  one  finds: 
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with  Xk  = X^k 
Aa(A^)  = Xkw\Xk) 


fk 


w'(- — ) 

A 


Ao(Aiv)  = Aivw'(Aiv) 


1 


w'(-p=) 


£ Ao(Xt)  = Xw'(X) 
k=0 


(«) 


Since  the  left  hand  side  of  the  equation  can  be  obtained  from  data  in  simple  extension  and 
compression  experiments,  the  function  w ’(A)  is  readily  obtained  and  equation  6 gives  a solution 
to  the  separation  problem. 

Here,  because  we  have  only  limited  compressive  data  and  substantial  extensional  data,  we 
chose  to  describe  the  stress-strain  response  of  the  rubber  with  the  semi-empirical  equation  of 
Martin,  Roth  and  Stiehler  (referred  to  as  the  MRS  equation  [10]),  which  is  known  to  fit 
adequately  data  of  cross-linked  polymers  [8,1 1]  for  stretches  up  to  X=3.5  and  for  compressions 
down  to  A=0.5  . The  MRS  equation  can  be  written  as: 

1 A (A-I) 

Ao(X)  = x (7) 

X 

The  two  adjustable  parameters  M and  A (nominally  material  constants)  can  be  found  by 
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fitting  the  tension  data  and  the  compressive  response  is  estimated  by  extrapolation  of  the 
function.  (Note,  in  future  work,  direct  compressive  experiments  will  be  performed  to  better 
determine  the  MRS  parameters.) 

The  Viscoelastic  Description  and  Bernstein,  Kearslv  and  Zapas  Theory 

In  1963,  Bernstein,  Kearsley  and  Zapas  [1,2,4]  proposed  a model  known  as  the  BKZ 
theory  which  describes  the  time-dependent  behavior  of  non-linear  viscoelastic  “perfect  elastic 
fluid’Mike  materials.  For  an  incompressible  material,  the  stress  ot]  (t)  can  be  calculated  for  any 
deformation  history,  using  the  BKZ  constitutive  equation: 


o.(r,7)  = -pb. 


f lU,(/l,/2,t-i,T)B:j(t,z)  - U2(I,J2,t-x,T)B:\t,x)}dx 
t 


with 


U = 

l 


dU 

V, 


(8) 


where: 
t is  the  time, 

T is  the  temperature, 
p is  the  hydrostatic  pressure, 

8 is  the  Kroneker  symbol, 

Bj/t,  t)  is  the  relative  left  Cauchy-Green  deformation  tensor, 
U is  the  strain  potential  function. 
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In  the  case  of  uniaxial  deformations,  the  relative,  left  Cauchy-Green  tensor  is  written  as: 


0 0 

fX(T)2  0 0 

\(x)2 

K& 

0 0 

, B~'(t,x)  = 

0 ^ 0 

W) 

X(t) 

/P 

o 

o 

0 0 ^ 

A(t) 

(9) 


and  the  deformation  invariants  are: 


MO2  + 2 Mr) 

a.(t)2  HO 


Krf  + 2 X(t ) 

\(t)2 


(10) 


A aft)  for  uniaxial  strain  is  obtained  by  the  equation: 


a o(o  = 2 

J 1 \{x)2 


X(t)  at/  + Kt)  du  dx 
k(f)  a/,  a(t)  a/2 


(ii) 


Single-step  deformation  history:  stress  relation 

Let  us  define  the  following  function: 
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(12) 


We  now  represent  the  single  step  stress  relaxation  in  simple  extension  : 


Ao  (f)  = 2 a2-I)(^;(r)  + |^(f)) 


(13) 


And  we  note  that  the  above  equation  resembles  the  Rivlin  and  Saunders  [12]  finite 
elasticity  equations,  except  that  Wn(t)  and  WI2(t)  are  time  dependent  as  are  I , and  I2.  Moreover, 
this  equation  applies  to  a material  which  is  not  assumed  to  be  in  its  equilibium  state,  since  the 
stress  may  be  relaxing.  That  is,  it  applies  to  viscoelastic  materials. 

Determination  of  responses  to  other  histories  from  single-step  stress  relaxation  data: 

Let  us  define  the  function  G(A,t,T)  as  the  value  of  Ao(t,T)  in  a single  step  stress  relaxation 
experiment  carried  out  at  a given  temperature: 


G(X,t,T)  =Ao(/\.,r,7) 


(14) 


and  its  partial  derivative  by: 


(15) 


We  can  then  predict  the  stress  history  in  any  uniaxial  deformation,  using  the  constitutive 
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equation: 


Ao(A,r,7)  = - f'G.(M,,-T,7)A 

J — A(l) 


(16) 


If  one  considers  deformations  for  which  the  material  is  at  rest  at  times  prior  to  x = 0 
(motion  following  the  rest  history),  equation  (16)  becomes: 


(17) 


Validity  of  the  time-dependent  BKZ  theory  model 

In  order  to  test  the  validity  of  the  VL-BKZ  approach  for  a given  material  and  stress  or 
deformation  histories  one  follows  the  sequence: 

1 . Perform  single  step  stress  relaxation  experiments  over  a range  of  deformations  and 
times,  at  a given  temperature. 

2.  From  the  relaxation  measurements,  determine  G,  (A,  t). 

3.  Perform  tests  in  arbitrary  deformation  or  stress  histories  and  compare  the  results 
calculated  from  the  BKZ  model  with  the  observed  responses; 


i.  In  multiple  step  stress  relaxation  or  constant  rate  of  deformation  histories,  one 
needs  to  simply  calculate  the  expected  response  from  equation  1 7 for  the  imposed 
deformation  history. 

i.i.  In  creep  experiments,  one  measures  the  deformation  and  inserts  the 
experimental  A(t)  into  the  right-hand  side  of  the  single  integral  equation,  and 
checks  if  the  calculated  stress  fits  the  experimental  stress.  For  example,  in  a 
constant  stress  experiment  does  one  calculate  a constant  stress  from  the  observed 
deformations? 
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III.  EXPERIMENTAL  METHODS[13a] 

Materials 


The  results  presented  in  the  subsequent  sections  were  obtained  for  a polyurethane  rubber 
designated  506  KW7  and  provided  by  the  NSWC,  Carderock  in  the  form  of  approximately  8 cm 
wide  and  2 m long  strips  of  4 mm  thickness.  The  material  was  not  otherwise  specified. 
Rectangular  samples  for  testing  in  uniaxial  extension  were  prepared  by  cutting  the  samples  with 
a guillotine-type  cutter.  The  samples  are  roughly  50  mm  long,  between  2 and  3 mm  wide  and  2 
mm  thick.  The  width  was  chosen  to  avoid  grip  effects.  That  is,  the  samples,  upon  testing,  were 
found  to  have  the  same  deformation  within  a 1 cm  gage  section  as  that  determined  by  measuring 
the  grip  displacement  and  dividing  by  the  original  grip  separation. 

Apparatus 

Environmental  testing  chamber 

In  order  to  validate  the  model  over  a range  of  experimental  conditions,  a PC  controlled 
humidity  and  temperature  chamber  was  designed  and  built.  It  is  depicted  schematically  in  Figure 
1 . The  temperature  is  controlled  by  an  integrated  circuit  which  receives  analog  signals  from  the 
computer.  The  temperature  probe  output  voltage  is  compared  to  the  computer  analog  signal.  The 
differential  voltage  is  amplified  by  a series  of  operational  amplifiers.  The  corresponding 
amplified  signal  triggers  on  or  off  a heating  device.  An  extra  cold  air  source  enables  the  chamber 
to  be  cooled  to  0 °C.  The  useful  temperature  range  for  the  apparatus  is  from  0 °C  to  80  °C  with  a 
measured  stability  of  ±0.2  °C. 
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In  a similar  manner,  a personal  computer  controls  the  relative  humidity  in  the  inner  box. 
Two  humidity  sensor  signals  give  a spatial  average  value  of  the  moisture  content  in  the  inner 
chamber.  The  signal  is  then  compared  to  that  which  is  proportional  to  the  desired  relative 
humidity.  Amplified,  the  differential  voltage  regulates  a proportional  integral  differential  (PID) 
electronics  system  which  activates  two  valves  to  provide  either  wet  or  dry  air  flow. 

The  relative  humidity  range  goes  from  5 % to  90  %,  with  a measured  stability  of  ± 2 %. 

Tension  and  compression  tests 

The  uniaxial  tension  and  compression  tests  were  performed  using  a load  frame  moved 
symmetrically  about  the  center  using  a stepper-motor  (Daedal,  Inc)  monitored  by  a personal 
computer  for  command  and  acquisition  of  input  signals.  The  system  is  described  in  Figure  2. 
Commands  are  sent  to  the  motor  through  an  indexer  and  a board  coupled  together.  Both  the 
velocity  and  the  acceleration  of  the  motor  screw  are  calculated  to  reach  the  desired  displacement 
or  stress  with  the  minimum  of  overshoot.  The  sample  length  is  measured  by  a miniature  Linear 
Variable  Transformer  (LVDT  type  MHR-100  or  MHR-500,  Schaevitz,  Inc.)  connected  to  an 
Analog  Transducer  Amplifier  (ATA  101,  Schaevitz,  Inc.).  The  signal  is  then  converted  into  a 
digital  one.  The  uncertainty  of  the  extension  ratio  is  less  than  0.1  %.  The  force  applied  is 
measured  using  a load  cell  (MB- 10  or  MB- 100  type,  Interface,  Inc.).  The  range  of  forces  studied 
goes  from  5 lbs  (25  N)  to  250  lbs  (1250  N) , depending  on  the  load  cell  features. 

The  true  stress  in  the  deformed  sample  is  calculated  assuming  that  the  material  is 
incompressible.  The  PC  reads  both  the  specimen  length  and  the  applied  force  through  an  AD/DA 
converter  card  (RTI-815,  Analog  Devices,  Inc).  The  extensometer  (respectively  load  cell)  signal 


12 


is  used  as  a feedback  to  control  the  strain  (respectively  the  stress)  in  stress  relaxation 
(respectively  creep)  experiments. 

A program  was  written  in  a Turbo  Pascal  code  to  control  the  apparatus  in  order  to 
perform  single  or  multiple  step  stress  relaxation  or  creep  tests.  Constant  strain  or  stress  rate 
experiments  can  be  carried  out  as  well. 

IV.  RESULTS 

Single-step  Stress  Relaxation  Experiments 

Typical  results  for  a single  step  stress  relaxation  experiment  are  shown  in  Figure  3a. 
The  uncertainty  of  the  experimental  stress  data  was  determined  to  be  approximately  5 % as  a 
single  standard  deviation  by  replicate  measurements  at  1 0 °C  with  samples  taken  from  arbitrary 
positions  within  a single  strip  of  material  provided  to  NIST.  An  F-test  [13b]  showed  that  the 
standard  errors  obtained  at  three  different  strains  (20  %,  50  %,  80  %;  0.20,  0.50,  0.80)  at  10  °C 
are  not  significantly  different  within  a probability  of  95  %. 

In  order  to  compensate  the  time  that  the  motor  takes  to  reach  the  desired  strain,  a 
systematic  shift  of  half  of  the  loading  time  is  applied  to  the  experimental  data  [14]. 

Variable  Separation  Assumption 

We  assume  that  the  response  is  time-temperature  and  strain  separable  and  use  separable 
functions  to  describe  the  stress  history: 
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(18) 


o(X,t,T,rh ) = gMRS(\).f[t,T).h(rh. ) 

with', 
t : time, 

T : '.temperature, 
rh  \ relative  humidity. 

Relative  Humidity  Effect 

Both  mass  up-take  and  mass  loss  experiments  are  currently  being  carried  out  on  the 
polyurethane  to  study  the  absorption  and  desorption  of  water.  Further  experiments  performed  in 
temperature  and  relative  humidity  controlled  conditions  will  enable  us  to  analyze  the  impact  of 
the  moisture  content  on  the  material  response.  The  relative  humidity  parameter  will  then  be 
included  in  the  model  developed  to  describe  the  polyurethane  resin  behavior. 

Time-temperature  Dependence 

The  time-temperature  dependent  function  f(t,  T)  was  obtained  using  data  from  stress 
relaxation  experiments  of  3000  seconds  duration.  The  range  of  temperatures  considered  was 
from  10°C  to  65  °C.  The  reference  temperature  T0  was  taken  as  10  °C.  The  horizontal  shift 
factors  were  calculated  from  a WLF-type  equation  (Eq.  19),  while  the  vertical  shift  is  given  by 
bTo(T)  = T0  / T.  An  uncertainty  of  ±10  % was  allowed  on  bTo(T)  to  optimize  the  curve  fit  to  the 
master  curve. 


Log(ar)  = 


c^T-T,) 

c2HT-T0) 


with  Cj  =9.57±0.1 1 , c2  =(123.7±0.9)  °C,  T0  = \0°C 


(19) 
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The  uncertainties  on  the  parameters  c,  and  c2  represent  a single  standard  error  of  estimate  from 
the  curve  fitting  procedure. 

The  master  curve  equation  for  the  stress  relaxation  response  was  then  fitted  to  a modified 
stretched  exponential  function: 


(20) 


Figure  3b  depicts  the  experimental  master  curve  and  predictions  from  the  fits  to  equation 
20  extrapolated  to  20  years  (7xl08  s).  To  give  an  idea  of  the  possible  uncertainties  in  such 
extrapolations,  we  have  chosen  to  also  present  the  bounds  assuming  that  there  is  a ±10  % 
variability  in  the  magnitude  of  the  modulus  (i.e.,  gMRS(A,)).  This  represents  approximately  2 
standard  deviations  on  the  mean  of  such  measurements  (see  first  paragraph  under  Results).  We 
ascribe  no  obvious  physical  significance  to  the  two  parameters  tc  and  (3  [15].  Furthermore,  we 
note  that,  to  our  knowledge,  a value  of  the  stretching  exponent  (~P)  that  is  less  than  0 has  not 
been  used  previously  to  describe  stress  relaxation  data. 

For  infinite  times,  the  right-hand  side  of  equation  20  approaches  1 and  captures  the 
network  behavior: 


lim  y?fr)  = 1 , with  tr  = — 


(21) 


The  parameters  tc  and  (3  are  determined  by  fitting  the  curve  using  the  Fletcher  and  Reeves 
Conjugate  Gradient  method  [16]  and  are  depicted  in  Table  1. 
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Table  1.  Curve  Fitting  Parameters  for  Equation  20.  (See  Text  for  Discussion). 


parameter 

value 

standard  error 

T 

2.14  105  s 

- 

P 

0.026 

5%  ! 

Moreover,  single  step  stress  relaxation  experiments  performed  on  the  506  KW7  elastomer 
at  1 0 °C  and  for  strains  up  to  1 00  % showed  that  the  time-temperature  function  f(t,  T)  does  not 
depend  on  the  stretch  ratio  X (for  1 <X<2). 

Figure  4 depicts  the  stress  relaxation  data  recorded  over  a time  scale  of  4 decades  along 
with  the  fit  using  the  modified  stretched  exponential  function  (Eq.  19).  The  experiment  was 
carried  out  at  10  °C.  For  times  greater  than  10  seconds,  the  function  f(t,T)  captures  reasonably 
well  the  time  dependent  process. 

Strain  Dependence: 

The  semi-empirical  Martin,  Roth  and  Stiehler  (MRS)  equation  was  used  to  fit  the  stress- 
strain  behavior  of  the  rubber.  The  MRS  equation  is: 

g^a)=Ao(A)=M(l-I)/a_I>  (22) 

To  determine  the  parameters  M and  A,  uniaxial  deformation  experiments  were  performed 
at  10  °C,  using  the  resin  type  D506  KW7.  Tension  results  are  depicted  in  Figure  5 as  the 
normalized  true  stress  versus  the  applied  stretch  ratio.  The  true  stress  was  normalized  by  dividing 
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the  experimental  values  by  the  reduced  time  function  (Eq.  20) f(tr,T).  Hence,  at  each  A,  three 
stress  points  are  plotted,  one  for  each  isochronal  time  value  (10  s,  100  s,  1000  s).  Also  shown  in 
Figure  5 is  the  least  squares  fit  of  the  MRS  function  to  the  isochronal  data  (10  s,  100  s,  1000  s). 
Note  that  the  fit  seems  to  deviate  somewhat  from  the  data  at  the  highest  stretches.  This  error  may 
be  due  to  a breakdown  of  the  time-strain  separability  assumption  or  to  the  time  (resp.  strain) 
dependent  power  law  (resp.  MRS)  function.  The  results  of  fitting  the  data  to  equation  22  are 
presented  in  Table  2.  Further  work  is  to  be  done  to  examine  other  possible  material  models  than 
the  MRS  model.  However,  for  the  current  purposes  the  MRS  function  is  adequate  to  proceed  to 
calculate  the  material  response  in  other  deformation  histories.  (As  an  aside  we  add  that 
preliminary  experiments  in  the  compression  mode  (0.8<A<1)  were  conducted  on  cylindrical 
samples  of  the  type  506  KW7  polyurethane,  but  the  results  could  not  be  used  due  to  technical 
problems.  A new  compression  apparatus  has  been  built  and  this  will  allow  us  to  perform  tests  in 
the  compression  region  down  to  a stretch  ratio  of  approximately  0.75.  Hence,  the  results  reported 
here  include  only  MRS  estimates  determined  from  tensile  data  alone.}  Also,  we  remind  the 
reader  that  the  purpose  of  using  the  MRS  equation  is  that  a means  to  extrapolate  into 
compression  is  required  to  predict  complicated  strain  histories  because  the  equations  in  the 
hybrid  VL-BKZ-type  model  require  relative  deformations,  which  become  compressive  when  the 
deformation  becomes  non-monotone  and  for  the  lateral  strains  in,  e.g.,  uniaxial  deformations. 
Table  2.  Results  from  Fitting  Tensile  Data  to  Equation  22.  (See  Text  for  Discussion). 


parameter 

value 

standard  error 

M 

4.33  MPa 

9% 

A 

0.07 

2% 
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The  parameter^  in  equation  22  was  found  to  be  +0.07  which  differs  somewhat  from  the 
range  0.3  - 0.4  reported  in  the  literature  for  cross-linked  rubbers  [17].  In  addition  to  the  visual 
quality  of  fit  observed  in  Figure  5,  one  can  get  a quantitative  estimate  by  examination  of  the 
normalized  residual  plot  shown  in  Figure  6.  While  the  maximum  relative  deviation  between  the 
fitted  curve  and  the  data  can  be  estimated  to  be  approximately  ± 6 %,  although  there  is  some 
structure  to  the  residuals,  there  is  little  systematic  deviation  for  the  different  times. 

Valanis  Landel  Function  Calculations 

Having  established  the  time  and  strain  dependent  functions  in  the  MRS  equation,  the 
calculation  of  w ’(A,t, T)  is  now  possible.  For  a given  stretch  ratio  X0  ,w’(Ag,t,T)  is  given  by: 

x]-\ 

x0wXkQit,D  ^ Vi  A(~r) 

-AO  k= 0 K (23) 

with  A,a.=A,0  !/2) 

A Turbo  Pascal  program  was  written  to  make  this  calculation.  In  practice,  only  10  terms 
are  necessary  to  obtain  convergence. 

VL-BKZ  Function  Calculations 

One  can  now  determine  the  Kernel  function  G,(A,t,T)  used  to  calculate  the  BKZ  single 
integral  for  a uniaxial  motion  following  the  rest  history. 
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(25) 


r n fT,  dG(X,t,T)  d 

with  GXkt,T)  = = — 

dt  dt 


1 /,  1 


Xw  (X,t,T)~ — w'{ — ,t,T) 
sfl  six 


The  Kernel  function  becomes: 


G (X  t T)-X  1 dw'{\lslX,t,T) 

dt  yx  dt 


with 


dw'(X0,t,T ) i 


dt 


X n 


k=0 


(26) 


d/m 

dt 


This  gives  the  stress  history,  using  the  single  integral: 


Ao(A  ,t,T)  = G(X(r),t,T)  - j 

0 


x(t)  dw'(X{tyx(x),t-x,T) 


X{x) 


dt 


\ 


X(x)  dwjy/X(x)/X(t),t-x,T) 


X(t) 


dt 


dx 


(27) 


Prediction  of  The  Stress  History  For  a Multiple-step  Uniaxial  Motion  Following  The  Rest 

History 

Using  the  experimentally  determined  single-step  stress  relaxation  response  function 
(equation  1 8),  one  can  predict  the  stress  (respectively  strain)  response  to  multiple  step  strain 
(respectively  stress)  histories.  These  are  examined  in  the  next  paragraphs. 


Stress  Relaxation  Experiment 

Multiple  step-up  stress  relaxation  experiment 

A triple  step-up  strain  history  was  applied  to  the  rubber  sample  (D  506  KW7)  at  a 
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constant  temperature  of  50  °C.  The  strain  and  the  true  stress  signals  were  recorded  during  3000  s 
after  each  deformation.  The  strain  history  (Figure  7)  is  summarized  by: 


A(x)  = 1 t<0 

A(x)  = Aj  0<x<tj 

A(x)  = A2  rj<x<T2 

A(x)  = A3  t2<x 

here 

Aj  =1.3  A2  = 1.4  A3  = 1 .5 
^=3000  5 r2=6000  s 


(28) 


The  BKZ  single  integral  (Eq.  27)  can  be  explicitly  calculated  and  the  triple  step-up  stress 
relaxation  response  is  given  by: 


0 <t<ty:  Ao(t,T)  = G(Aj,f,7) 


0<'<V 


Ao(t,T)  - G 


(k  1 

' A, 

— t-t  T 

+ G(X2,t,T)  - G 

— ,t,T 

J 

A. 

V 1 / 

(29) 


[ A3  1 

[ ^3  l 

( A,  ] 

U ) 

t2<t:  Ao(t,T)  = G 

— t-t  T 
A ’ 2’ 

\ 2 / 

+ G 

— t-t  T 

a/  ” 

\ 1 y 

- G 

— ,t-t.,T 

A,  ! 

\ 2 y 

+ G(  A3 , 7)  - G 

— t T 

A, 

V 1 / 

We  remark  that  Equation  29  applies  to  ideal  steps.  In  Figure  8 we  depict  the  material 
response  to  the  three  step  deformation  history  along  with  the  BKZ  computation  based  on 
Equation  29.  The  agreement  between  model  and  experiment  shown  in  Figure  8 is  reasonable.  In 
Figures  9-1 1 a comparison  between  the  calculated  values  and  the  measured  values  is  shown  for 
the  individual  steps  with  the  time  reset  to  the  beginning  of  each  step.  These  figures  show  that  the 
BKZ  calculation  from  Equation  29  is  reasonable,  though  the  calculated  short  time  response 
deviates  systematically  from  the  measured  response.  The  first  reason  for  this  is  the  finite  rate  of 
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deformation  effect  that  is  seen  in  the  measured  responses.  As  mentioned  previously,  this  could 
be  corrected  by  taking  away  half  of  the  loading  time  (several  seconds)  from  the  experimental 
results  [14].  This  is  not  done  here  since  the  purpose  is  to  illustrate  that  the  simple  BKZ-step 
calculation  provides  a reasonable  estimate  of  the  multiple  step  history.  In  all  subsequent 
calculations  we  integrate  numerically  the  full  strain  history  (Equation  27)  rather  than  use 
Equation  29. 

Multiple  cycle  stress  relaxation  experiment: 

The  506  KW7  material  stress  response  to  a quadruple  cycle  strain  history  (Figure  12)  is 
shown  in  Figure  13  along  with  the  hybrid  VL-BKZ  simulation.  The  experiment  was  performed  at 
10  °C.  The  stress  relaxation  computation  was  made  using  Equation  27.  While  the  calculation 
does  not  fit  the  stress  response  for  the  unloading  steps,  it  does  describe  reasonably  well  the 
relaxation  process  corresponding  to  the  loading  step  of  the  four  cycles.  This  is  as  anticipated 
because  the  BKZ-type  constitutive  models  are  known  to  provide  less  good  predictions  when  the 
sample  is  unloaded  or  the  load  is  reversed  [3], 

A point  to  be  made  here  is  that  the  ‘unload’  portions  of  the  deformation  history  were  not 
fully  unloaded.  This  is  because  a full  return  of  the  sample  to  a nominal  zero  deformation  caused 
sample  buckling  (due  to  compression).  Therefore,  the  deformation  was  cycled  between  the 
maximum  extension  and  a small  positive  stress  to  avoid  the  buckling  phenomenon. 

Creep  Experiments 
Multiple  step-up  stress  history 
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A creep  experiment  was  performed  by  applying  a four  step-up  stress  history,  at  a 
temperature  of  10  °C.  The  stress  levels  were  chosen  to  cover  the  deformation  range  of  >2=1.2  to 
>2=1.7  as  seen  in  Figure  14.  Each  step  lasted  2500  seconds.  The  measured  strain  history  A(t) 
observed  during  the  multiple  creep  steps  and  shown  in  Figure  14  was  used  in  Equation  27  to 
compute  the  stress  history.  Figure  15  compares  the  applied  stress  history  with  that  calculated 
from  the  hybrid  VL-BKZ  model.  The  model  gives  a relatively  good  fit  to  the  experimental  data 
corresponding  to  the  four  stress  steps  (Figures  15  through  19).  We  do  note  that,  while  the 
deviations  between  the  model  calculations  and  the  data  are  less  than  about  1 0 % of  the  applied 
stress,  the  model  exhibits  a somewhat  softer  response  that  approaches  the  experimental  value  at 
long  times. 

Trapezoidal  stress  history: 

Figures  20  (respectively  21)  show  the  material  response  (resp.  the  deformation)  to  a 2 
trapezoidal-shaped  step  stress  history.  The  experiment  was  performed  at  23  °C,  using  a sample 
type  506  KW7.  The  model  is  within  reasonable  agreement  (less  than  10  % difference)  with  the 
experimental  data  recorded  during  the  increasing  stress  ramps.  As  observed  for  the  four  cycle 
stress  relaxation  experiment,  the  model  does  not  describe  the  unloading  sections  of  the  history  as 
well  as  it  does  the  loading  ones. 

V.  CONCLUSIONS  AND  FURTHER  WORK 

A time  dependent  Valanis  - Landel  strain  energy  function  has  been  used  to  describe  the 
non-linear  response  of  a polyurethane  rubber.  The  time  dependent  VL  function  was  determined 
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by  performing  uniaxial  single  step  stress  relaxation  experiments,  assuming  that  the 
incompressible  material  response  is  time  and  strain  separable.  Experiments  were  performed  over 
a temperature  range  of  5 °C  to  70  °C  for  strains  from  5 % to  100  %.  The  time  - temperature 
superposition  function  was  established  using  data  obtained  at  X=\  .5  and  for  temperatures  ranging 
from  10  °C  to  65  °C.  The  time-temperature  dependent  strain  energy  function  was  then  used  in  the 
framework  of  the  BKZ  single  integral  viscoelastic  constitutive  equation.  The  model  developed 
predicts  the  non-linear  viscoelastic  material  response  to  multiple  cyclic  or  step-up  stress  or  strain, 
and  constant  rate  of  deformation  histories  in  uniaxial  extension  within  an  error  of  less  than  1 0 % 
for  the  range  of  parameters  examined. 

Ongoing  experiments  to  investigate  the  impact  of  moisture  on  the  material  response  will 
help  in  the  prediction  of  the  behavior  of  the  polyurethane  rubber  under  various  environmental 
conditions.  Initial  modeling  will  include  the  relative  humidity  variable  in  the  time-temperature- 
strain  function  for  the  rubber  as  a moisture  shift  factor,  similar  to  that  used  in  the  time- 
temperature  superposition. 

In  addition,  new  experiments  to  obtain  reliable  data  in  compression  will  enable  us  to 
describe  both  the  tension  and  compression  isochronal  stress-strain  results  with  a better 
agreement.  Other  models  than  the  Martin,  Roth  and  Stiehller  equation  will  be  investigated  to  fit 
the  stress-strain  isochrone  curves,  such  as  the  Gaylord-Douglas  model  [18].  Furthermore,  other 
constitutive  equations  describing  solid-like  material  responses  found  in  the  literature  will  be 
compared  to  the  Bernstein,  Kearsley  and  Zapas  model. 
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Figure  1.  Schematic  of  Environmental  Testing  Chamber  (See  text). 
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Figure  2.  Schematic  of  Control  System  for  Stepper  Motor.  (See  text.) 
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Figure  1.  Schematic  of  Environmental  Testing  Chamber  (See  text). 
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Figure  2.  Schematic  of  Control  System  for  Stepper  Motor.  (See  text.) 


Figure  3a:  Typical  strain  (left)  and  stress  (right)  versus  time  curves  for  a 
single  step  stress  relaxation  test  on  a polyurethane  elastomer. 


Peak  strain  = 50%  ; T = 10°C. 
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Figure  3b:  time  - temperature  superposition  master  curve 

X=  1.5. 
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Figure  4:  Stress  versus  Log(time)  in  stress  relaxation  test  at  X=1.5  and  T=10°C. 
The  modified  stretched  exponential  is  given  by  Eq.  (19). 
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Figure  5:  Comparison  of  normalized  stress  function  vs  deformation 
in  tension  with  curve  fit  response  given  for  the  MRS  equation. 

T =10°C. 
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Figure  6:  Residuals  between  normalized  stress  function  (see  text) 
and  the  MRS  curve  fit.  Data  are  from  Figure  5. 
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Figure  7:  Strain  versus  time  curve  for  a multiple  step-up  stress  relaxation  experiment. 
T=50°C. 
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Figure  8:  Stress  versus  time  curve  for  a multiple  step-up  stress  relaxation  experiment 
T=50°C. 
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Figure  9:  Log(stress)  vs  Log(time)  curve  for  a multiple  step  stress  relaxation  experiment. 
T=50°C. 
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Figure  10:  Log(stress)  vs  Log(time)  curve  for  a multiple  step-up  stress  relaxation  experiment. 
T=50°C. 
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Figure  1 1 : Log(stress)  vs  Log(time)  curve  for  a multiple  step-up  stress  relaxation  experiment. 
T=50°C. 
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Figure  12:  Strain  versus  time  curve  for 
a multiple  stress  relaxation  cyle  history. 
T=10°C. 
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Figure  13:  Stress  versus  time  curve  for  the  material  response 
to  multiple  stress  relaxation  cycles. 

T=10°C. 
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Figure  14:  Strain  versus  time  curve  for  a multiple  step-up  creep  experiment. 
T=10°C. 
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Figure  15:  Stress  versus  time  curve  for  a multiple  step-up  creep  experiment 
T=10°C. 
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Figure  16:  Log(stress)  vs  Log(time)  curve  for  a multiple  step-up  creep  experiment. 
T=10°C. 
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Figure  17:  Log(stress)  vs  Log(time)  curve  fora  multiple  step-up  creep  experiment. 
T=10°C. 
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Figure  18:  Log(stress)  vs  Log(time)  curve  for  a multiple  step-up  creep  experiment. 
T=10°C. 

STEP  3. 
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Figure  19:  Log(stress)  vs  Log(time)  curve  for  a multiple  step-up  creep  experiment 
T=10°C. 

STEP  4. 


Figure  20:  strain  versus  time  vurve  for  a material  response 
to  a trapezoidal-shaped  stress  history. 
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Figure  21 : stress  versus  time  curve  for  a material  response 
to  a trapezoidal-shaped  stress  history. 

T=23°C. 
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